A simple solvable energy landscape model that shows a thermodynamic phase 

transition and a glass transition 
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When a liquid melt is cooled, a glass or phase transition can be obtained depending on the cooling 
rate. Yet, this behavior has not been clearly captured in energy landscape models. Here a model is 
provided in which two key ingredients are considered based in the landscape, metastable states and 
their multiplicity. Metastable states are considered as in two level system models. However, their 
multiplicity and topology allows a phase transition in the thermodynamic limit, while a transition 
to the glass is obtained for fast cooling. By solving the corresponding master equation, the minimal 
speed of cooling required to produce the glass is obtained as a function of the distribution of 
metastable and stable states. This allows to understand cooling trends due to rigidity considerations 
in chalcogenide glasses. 

PACS numbers: ***** 



Humankind has been using glassy materials since the 
dawn of civilization. However, their process of formation 
still poses many questions _1][2][3][4] 5] 6] 7]. Glasses do 
not have long range order and are formed out of thermal 
equilibrium, resulting in a limited use of the traditional 
tools of the trade in solid state and statistical mechanics. 
Moreover, numerical simulations are not able to provide 
definitive answers, since cooling speeds achieved in nu- 
merical simulations are orders of magnitude higher than 
in real cases [8]. One of the main issues is the nature 
of the glass transition [9|], for example, is it a purely dy- 
namical effect or there is a underlying thermodynami- 
cal singularity? The answer to this question has practi- 
cal implications, as how to calculate the minimal cool- 
ing speed depending on the chemical composition in or- 
der to form a glass, or why some chemical compounds 
form glasses while others will never reach such state [lO| • 
Concerning this relationship between chemical compo- 
sition and minimal cooling speed, Phillips [lO| observed 
that for several chalcogenides, this minimal speed is a 
function of rigidity. This initial observation was the ing- 
nition spark for the extensive investigation on rigidity of 
glasses [nj [12 [11 [13 [H [3 [13 El, yet this basic obser- 
vation has not been quantitatively explained. 

On the other hand, the energy landscape has been a 
useful picture to understand glass transition '^. How- 
ever, due to its complicate high dimensional topology, 
it is difficult to obtain closed analytical results. It is 
not even clear how a phase transition is related with the 
topology of the landscape, i.e., why a global minimum 
leads to singularities in the thermodynamical behavior. 
Clearly, there is a lack of a minimal simple solvable model 
of landscape that can display a phase and a glass transi- 
tion depending on the cooling rate. Here we present such 
model by combining the two most basic ingredients that 
are belived to be fundamental in the problem. Further- 
more, the model allows to get a glimpse on the connec- 
tion between minimal cooling speed, energy landscapes, 
rigidity, and Boolchand intermediate phases |19||2(j|. 



The first ingredient is based in a well known fact: 
glasses are trapped in metastable states, while crystals 
are global minimums in the landscape. A common way 
to describe the corresponding physics is through the use 
of two level system (TLS). If the glassy metastable state 
has energy Ei and the crystalline global minimum an en- 
ergy i^o — 0, the system is trapped in the glassy state 
due to an energy barrier V measured from Ei, as seen 
in Fig- JE Following Huse et. al.[21| and Langer et. 
al. 122] [23] [24'| , who described the residual population of 
the metastable state for a TLS at zero temperature, the 
cooling process can be described by a master equation, in 
which the probabibilty p[t) of finding the system in the 
metastable state, assuming that the system is in contact 
with a bath at temperature T, is[22|. 



dp{t) 
dt 



-rt;pW + r;t[i-pW] 



(1) 



where T^^ is the transition rate from the upper well to 
the lower, and the transitions from the lower to the upper 
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take place at rate T ^-^ = e t T^^. If quantum mechan- 
ical tunneling is neglected, V^^ — Foe^", where Fq is a 
small frequency of oscillation at the bottom of the walls. 
Eq. (H)) describes the relaxation towards po{T), the 
population at thermal equilibrium obtained from the sta- 
tionary condition, as can be seen by rewriting Eq. ^ as 



dpjt) 
dt 



^Tn{l + e-^){po{T)~p{t)), 



where Po{T) is given by, 

Po{T) = 



e T 
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(3) 



When the system is cooled by a given protocol T = T{t), 
it can be proved that at zero temperature there is a prob- 
ability p{T = 0) for the system to be in the metastable 
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FIG. 1: The two level system energy landscape, showing the 
barrier height V and the asymmetry Ei between the two lev- 
els. The population of the upper well is p{t). 



state, which is indicative of a glassy behavior |23l|[25l|. 
This simple model is very appealing and can be used 
to explain low temperature anomalies in glasses [2 6] 27]. 
However, a huge part of the physics is missing: the sys- 
tem does not present a phase transition at low cooling 
speeds. To achive this goal, here we introduce a key el- 
ement to the TLS landscape topology: the multiplicity 
of states. Again, there is a common agreement that the 
number of metastable states is much bigger than their 
crystalline counteparts. Assume that the energy Ei has 
a degeneracy gi, while the ground energy Eq has degen- 
eracy go, thus Eq. ([T]) needs to be modified to take into 
account transitions between different states that are in 
the low and upper wells. Call p^s(i) the population of 
one of these gi in the upper states, and Pisit) the pop- 
ulation of one of these go in the low states. Eq. ([T]) 
becomes. 
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So 



+ E A^Mt) + E rr/Pi™(^). (5) 

where TJi denotes the transition rate from state s to Z, 
both in lihe upper well. The notation for the other tran- 
sition rates is similar, and an equivalent expression can 
be written for dp^s{t)/dt. To formulate the model, we 
use the simplest topology, i.e., all metstable states are 
connected within them with the same transition rate, 



I.e. , -L -K^- 

talline states Ff', 



^-|-|-. A similar situation holds for the crys- 
^4,4,. Transitions between up and 



lower states have also the same probability F; 



and Fl',, 
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-n 



-it- 



Under such simple landscape topology. 



the previous master equation can be reduced to. 



dpjt) 
dt 



-goT^nPi^) + 5ir;t [i - pW] 



(6) 



where now p{t) = ^ Ptsi^) is the total probability of 
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finding the system in the upper well. Let us show how Eq. 
([^ can give a phase transition under thermal equilibrium 
conditions. In that case, dp{t)/dt = and, 



Po{T) = 



.gir^t 



^ (gi/go)e ^' 



(7) 



A phase transition can occur if {gi/go)e^^^ becomes 
discontinous in the thermodinamical limit. The most 
simple example is the following. Suppose that we have 
N particles, and the potential is such that the crystalline 
state is unique {go = 1), with energy Eq = 0, and assume 
that the number of metastable states grows exponentially 
with N, as is the case in many glassy systems Q where 
gi = e^'^^^^^i). 51(£'i) is a measure of the landscape 
complexity [9|. Also, the only way to make (E) an inten- 
sive quantity with only one energy is to have Ei — Ne, 
where e is an energy per particle. As an example, this 
behavior can be readly obtained when two particles, con- 
fined in cells, interact with neighboring cells as in nearly 
one dimensional models of magnetic walls [2^ . For this 
particular case, gi — 2^ and go = 1. Using the previous 
general considerations, po{T) can be written as. 



PoiT) 



^N 



^{lnniEi)-j,)N 



_^g(lnO(Bi)-f)jV 1 + z 



N 



(8) 



with z = [cxp(n(Ei) — e/T)N]. In the thermodynamic 
limit N — )■ oo, the function f{z) = z^ develops a dis- 
continuty at z — 1, and it is easy to see that there is a 
phase transition at temperature. 



T,^ 



lnn{Ei) 
with a discontinuous specific heat, 

dpoiT) _\ OifT^Tc 



c = e- 



dT 



HT ^Tr 



(9) 



(10) 



Now the model is able to produce a phase transition un- 
der thermal equilibrium. This can be clearly seen in Fig. 
2 for gi — 2^ and go = 1, where we plot Eq. ([5]) for 
different values of N using dotted lines. Notice how the 
phase transition is built by a progressive sharpening of 
the jump in po{T) as N grows. According to Eq. PUI) . 
the specific heat is just the derivative of po(r), thus the 
sharpening leads to the singularity in the thermodynam- 
ical limit. 

Now we will show that a glassy behavior is obtained 
for fast enough cooling. To solve Eq. ([5]), one needs 
to specify the cooling protocol T = T{t), and write the 
master equation in terms of a dimensionless cooling rate. 
Two kinds of protocols are useful 22] [2^ , one is the linear 
cooling T = To — rt, used mainly in experiments, and the 
hyperbolic one T = Tb/(1 -I- Rt), which allows a simple 



handling of the asymptotics involved. For the hyperbolic 
case, the master equation can be written as, 



.dp{x) 
dx 



-gix^" + {go + gix'')p{x) 



(11) 



where x = exp(-V^/r) and 5 = RV/ToTq. The parame- 
ter fx = Ei/V measures the asymmetry of the well. The 
linear case also follows Eq. PT|) . since one can rescale the 
boundary layer [2^ that appears in Eq. PTj) . leading to 
the same hyperbolic equation with 5 = rV/T^. Eq. pTJ) 
can be solved to give, 



p{x) = exp 

|p(o) 
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(12) 



As an example. Fig. 2 shows p{x) for different cooling 
rates and system sizes, using a linear cooling and gi = 
2^, go = 1, compared with the equilibrium distribution 
that develops a phase transition at Tc- Notice in Fig. 2 
that p{0) is the residual population at T = 0, indicative 
of a glassy behavior. Also, the slope of dp{T)/dT does 
not tend to infinity, and the corresponding specific heat 
c in no longer discontinuous, as in real glass transitions. 
We can obtain the analytical value of p(0) by assum- 
ing that the system was at thermal equilibrium before 
being cooled at a temperature Tq » Tc- In that case 
a; — >■ 1, and the population is given by the equilibrium 
distribution, poixo) = {gi / go)x'^ / {{gi / go)xo + 1) where 
xo = exp(— T^/Tq) ■ From Eq. (fTTj) . we obtain a general 
expression for p(0). 



p(0) = 



2^0 + (50/51) 
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Zero population is only achieved if both terms in Eq. 
(fT5|) are zero, as is the case for 5 — >■ 0. Then we recover 
the phase transition, a fact that makes us confident in 
the result. To understand more deeply Eq. (TTSl) . let us 
study the particular case gi = e^'"^(^i) and go = 1, 
with El — Ne. The second integral contains the term 
g^y'^+^' k. exp [(lnrj(£'i) - e/T)N], which can be or 00 
in the thermodynamical limit depending whether y < Xc 
or y > Xc, where. 



exp(-F/T,) - n{Ei) 



^V/e 



(14) 



If V does not scale with N, for big N Eq. (IT51) can be 
written as, 



p(0) 



2^0 + (50/51) 



5i< 
6fi 



exp -^^Y-^ +gij{l + p,,Xc/S)6'^, 

(15) 

here 7 is the lower incomplete gamma function. The 
evolution of the residual population given by Eq. (J15p 




FIG. 2: Population as a function of the temperature using 
a linear cooling for different number of particles A'' = 2, 4 
and 8, with V = 1.0, e = 1, i? = 1.4 and Tq/Tc = 72 , 
gi = 2^ , go = I obtained by solving the master equation. 
The equilibrium population, obtained for 5 — >■ is also shown, 
as indicated in the inset. Notice how the phase transition is 
built by a progressive sharpening of the jump in p{T) as TV 
grows. 



is shown in Fig. ([3]) as a function of the cooling speed 
and system size. As a general trend, the cooling speed 
required to make a glass with fixed p(0) increases with 
the system size. Also, it is possible to observe a crossover 
which separates different behaviors of p(0). For example, 
in Fig. ([3]), if iV = 4, 5, 6 and 7, p{0) begins to increase 
for a high S after it reaches a plateau that begins around 
(5 « 1. The same increase is observed for iV = 9 and 10, 
although shifted to the right in such a way that it does 
not appear in the current plot. This crossover is due to 
the different growing speeds in Eq. dTSl) . The first term 
of Eq. pS]) goes to zero if. 
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SciN), 



(16a) 



which defines a critical speed Sc{N). For 6 > Sc{N), p(0) 
is dominated by the first term in Eq. (J15p . The remaining 
term in Eq. (jlSp regulates the residual population for 
lower speeds S < 6c{N). This term produces the plateau 
at a saturating value of p(0). 



p(0) 



1 



1 + (Ne/V) 



Ps{0). 



(17) 



For a finite N, this implies that there are two kinds of 
glassy phases, one obtained for intermediate cooling rates 
in which p(0) reaches a limiting value. The other kind is 
obtained for 6 > 6c{N). 

However, since Ps{0) goes to zero as iV ^> cxd, 6c{N) 
turns out to be the minimal speed required to make a 
glass in the thermodynamical limit, and leads to a critical 
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FIG. 3: Residual population at T = as a function of 
the cooling speed for different number of particles A'^ = 
10, 15, 20, 25 and 30, with V = 0.5, e = 1, i? = 1.4, ^i = 2^^, 
go = 1 and Tq/T, = 72. 
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(18) 



From an analysis of Eq. ([2]) and Eq. ([T]), it is easy to 
see that Re is basically the inverse relaxation time for 
crystalization. Although it is surprising that Re does 
not depend on V, this is a result of the assumption that 
V does not scale with N. If this is the case, the term 
giXQ ^ in Eq. (|13p can determine if goXo plays a role 
in the first exponential, leading to a critical speed that 
depends on V. Notice that the result is in agreement 
with the remarkable observation made by Phillips con- 
cerning chemical composition and minimal cooling speed 
required to make glasses [lOJ, since Re depends on the 
number of metastable states. Rigidity provides such an 
indirect count of metastable and stable states [2^ [30| . 

In conclusion, we have introduced the topology of the 
energy landscape in a two level model of glass. As a 
result, we have a solvable model that has a thermody- 
namic phase transition for low cooling rates and a glass 
transition for fast cooling. 
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